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Abstract 

We consider the transitions between electromagnetic bound states, such 
as the exclusive weak decay of a muonic atom into an electronic atom: 
(fj,~Z) — > (e~ Z)v e Va- We show that relativistic effects in the atomic wave- 
functions are crucial for determining the rate. In the case of heavy atoms, the 
exclusive channel branching ratios exceed 10 -6 , possibly bringing the study 
of these rare decays within experimental reach. Such processes thus provide a 
detailed laboratory for studying the high momentum tail of wavefunctions in 
atomic physics; in addition, they provide a simple toy model for investigating 
analogous exclusive heavy hadronic decays in quantum chromodynamics such 
as B — > irev. 
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1.) Introduction 

The calculation of the rate of a weak decay of a heavy hadron into an exclusive 
channel, such as B —>■ Dev or B — > jK*, poses a challenging problem in nonpertur- 
bative quantum chromodynamics because all the complexities of the hadronic wave- 
functions enter. In this paper we point out the possibility of studying an analogous, 
but far simpler process in atomic physics: the weak decay of one electromagnetic 
bound state into another. Simple examples occur in the decay of muonic atoms, 
such as 



Because the atom remains in a bound state, but the atomic or nuclear species is 
changed, we refer to such processes as "atomic alchemy." We could also consider the 
decay (n + e~) — > (/i + e~)z/ /t , or the decays of Coulomb-bound hadronic atoms such 
as (tt~ Z) — > (fi^Z)]?^ or — > (pZ)tt ; however, the observation of the latter 
transitions appears to be impractical because of the high rate at which hadrons are 
absorbed by the nucleus. 

Other interesting examples of atomic alchemy occur when a nucleus is forced 
to change its state because of external scattering, while an atomic electron remains 
bound. The scattering may be due to neutron-induced fission, nuclear Compton scat- 
tering, or photodisintegration; the latter is particularly interesting for deuterium, 
7((i + e~) — > (p + e~)n. One can also study atomic alchemy when a nucleus changes 
its state by an internal process, such as /5-decay. 

In this paper, we consider the exclusive transition B\ — > B 2 + X, where both B\ 
and B 2 are bound states consisting of particles (bis) and (62s), respectively, and the 
transition proceeds by the weak decay b x — > b 2 +X. The particle s will be called the 
spectator. We assume that there is an empty bound state in B 2 for the particle b 2 
to settle in; in the detailed calculations given below we take this to be the ground 
state. 

The basic features of atomic alchemy can be easily understood using non- 
relativistic mechanics; the analysis is similar to that for nuclear collisions involving 
capture of atomic electrons, as discussed, for example, by Migdal |]]. We review the 
major points before plunging into the full relativistic calculation. Because the weak 
decay occurs over a time short compared to the period of an orbit, we can use the 
"sudden approximation" . The probability amplitude for the atomic transition then 
factors into the free matrix element of the weakly decaying (moving) particle times 
a form factor F(q 2 ). This form factor is just the overlap of the wavefunctions of the 
initial and final states, which we write as ipi and ip 2 . In the rest frame of the initial 
state, F(q 2 ) is 



(/1 Z) -> (e Z)u e v, 
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Here k\ is the momentum of bi (see Fig. 1), and m re ^2 and v re i,2 are the reduced 
mass and the relative velocity of the final state particles b 2 and s. The velocity v re i j2 
is a function of q, the momentum transfer carried by X, which by conservation of 
momentum is equal to the recoil momentum of B 2 . We have 

kx-q h h q 

Vrel,2 = 1 = • (3) 

m b2 m s m red! 2 m b2 

In the decay (Z}i~) — ► (Ze~)v^v e) the argument k 2 of the wavefunction -02 is 
approximately k 2 ~ k\ — q. In momentum space the muon wave function does not 
vary much over the domain where the electron wave function differs appreciably 
from zero. The integral in Eq. (0) can therefore be approximated by 

F(g 2 )=^(^ 2 *(f=0) . (4) 

For small Z the momenta of both the muon and the electron are nonrelativistic and 
the effects of the finite nucleus size are small; we therefore can use for ipi and ip2 the 
Schrodinger wavefunctions for a point nucleus. For the IS states 

/ — 3/2 -1 

Uk)= r Vn \ l2 , a t = —,(1 = 1,2), (5) 
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where denote the Bohr-radii of the atoms. Because ipi(r = 0) = [^af^ 2 ^] x , 
Eq. (HI) becomes 



F(q 2 ) = 8 — 



\ 3/2 r ' 1 



m^J \l + q /(Zam 



A*/ J 



(6) 



The physics of capture into a hydrogenic nS state is controlled by the ratio of 
the change in nuclear velocity to the atomic velocity Za/n. Equation (|) suggests 
that the highest rates for the transition of muonic to electronic atoms will be found 
for heavy nuclei with [Zam^ 2 > q 2 . For these atoms, however, Eq. (^]) provides 
only a very rough estimate for the form factor because relativistic effects (mainly 
affecting the wavefunction of the {Ze~) atom) become crucial and greatly enhance 
the transition rate. 

In the case of the transition (7r + e~) — > (fi + e~)v fJ- the monoenergetic weak decay 
changes the species and velocity of the heavy particle serving as a nucleus. Since 
the momentum transfer from the pion to the muon is much larger than either Bohr 
momenta, the integrand in the form factor in Eq. (|2|) is dominated by those values of 
k\ for which either k\ or k 2 matches a Bohr momentum. Both contributions turn out 
to be equally important, and using for the IS states the Schrodinger wavefunctions 
for a Coulomb potential, we obtain for the form factor 

F(q 2 )~ (r = 0)ip* 2 (m red , 2 v rel , 2 )~ 16 f — ^ . (7) 

V V re l,2 J 
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Here f re z,2 = (m^ — m^)/(2m ir m )1 ) = 0.28, which is small enough that our non- 
relativistic treatment is justified. The probability for capture of the electron from 
the muonic decay is unfortunately very small: 

P = [F(q 2 )} 2 ~ 25Qa 8 /vf, el2 ~ 5 x lO" 11 . (8) 



2.) Relativistic Analysis of Atomic Transitions 



We now turn to a more detailed analysis which allows us to treat decays such 
as (Zjj,~) — > (Ze~)v fl i> e , where for large Z the velocities of the bound constituents 
are relativistic. In principle, the analysis requires a fully covariant description, 
such as the Bethe-Salpeter equation or the light-cone Fock state expansion. In 
Ref. a covariant description was given, but at the expense of a somewhat artificial 
momentum-dependent mass for the decaying constituent. While that description 
works well when the decaying constituent is much heavier than the spectator, it 
fails when the spectator is the heavier particle, because the high momenta of the 
constituents, which prove crucial, are cut off (see Eq. (2.1) of Ref. @). Here we 
use a non-covariant description that is adequate for the cases of interest, where 
the spectators in the initial and final bound states have small relative velocity, 
but where the velocities of the particles bound to the spectators can be relativistic. 
This description is particularly appropriate for the transition (/jTZ) — > (e~Z) +h , tl v e 
because the mass of a nucleus is much greater than m^. 

We begin with a relativistic treatment of the simple, monoenergetic, two-body 
atomic transition (7r + e~) — > (fi + e~) + i> fl . This is a close analog to the hadronic decay 
B — > Dev. The relative velocity of the (7re) and (/ie) systems is only v re \ = 0.28, 
so one can simplify the kinematics by assuming them to be at relative rest. We 
consider only S-wave bound states, so the spin of the bound state B\ is just the spin 
of the electron, and the final bound state B2 is either a pseudoscalar or a vector 
particle with the corresponding well-known spin combinations of the constituents. 
For example, the initial state B\ with spin projection R is represented in its rest 
frame as 

\B l ,p B = ^R) = J2m Bl / Mki)a+(h)b + (k 3 )\0) , (9) 

J ■ >^k\kl (2tt) 3 



where = —ki, and k® = (kf + m 2 ) 1 ^ 2 for (i = 1,3), and a + (ki) and b + (k 3 ) are 
creation operators for the constituents that act on the vacuum state |0). The state 
is "covariantly" normalized in the volume V, so that 

/d^ ki — * 
—J^^oi^i . (10) 
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The matrix element for the decay (7r + e ) — > (fi + e ) + is then written as Q 

AG F V ud r f d 3 h tptifa) ip%(k 2 ) a Um n 

M rs = n=— \ 4m Bl m B2 / i , S rs — — (11) 

V2 V J (2tt) 3 y^o 2 

5 rs = u r (g) 7 o^^^(fci-g) , fc 2 = K + (fci-<f) 2 ) 1/2 , (12) 

where /„- ~ 130 MeV is the pion decay constant, and V u d is the relevant CKM ma- 
trix element. Here u r (q) and v s (ki — q) are the Dirac spinors of respectively the 
neutrino with spin r and the muon with spin s. The argument k 2 of the bound 
state wavefunction ^ is given by k 2 = k\ — {m re d, 2 lm^)q. Note that in S rs we have 
kept only the zeroth component of the weak current since the pion is essentially at 
rest. In the limit k\ — <f <C m^, only the large components of the spinor need to be 
retained, and S rs takes the simple form 



e = ia 2 , = (1,0), x + = (0,l). (13) 

We square the matrix element, sum over pseudoscalar and vector final states, and 
average over the spin of the initial state. Using 

J2SrsS; s = 2m,\q\ , (14) 

rs 

the spin-averaged matrix element becomes 



\M\l = AGl\V ud \ 2 flml\F{q 2 )\ 2 m,\q\ 
F(q 2 ) = J j^A(ki)r 2 02) , (15) 



where we have used (see Eq. ([11])) the approximation Ak\k 2 ~ 4mB 1 m B2 . The 
two-body kinematics fixes the momentum transfer \q\ to be 

W\ = (™> 2 Bl -m 2 B2 )/(2m Bl ) « (ml -mj)/(2m 7r ) . (16) 

The decay rate in this approximation is 

2 2 

ml — m, 



To get an estimate of the errors made by using non-relativistic mechanics, we 
calculated the decay rate of a free pion making the same approximations, i.e. just 



2 We use spinors normalized as u r (p) u s (p) = —v r (p) v s (p) = 2mS rs 
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taking the zeroth component of the weak current and retaining only the large com- 
ponents of the Dirac spinor of the muon. We find 

rapP ro X _ G 2 F \V ud \ 2 flm^m^ml - ml) 2 

~ Mm% • (18) 

This differs from the rate calculated using relativistic mechanics only by the factor 

rjg? /YfZ° x = m^/m, « 1.3 , (19) 

which indicates a possible error of about 30% . The branching ratio for decay to the 



IS state, obtained by dividing Eq. ( O) by the approximate form in Eq. (18 



is 



BR((7r + e -)^(^ e -)^) = |F(g 2 )| 2 . (20) 

Using the Schrodinger wavefunctions of Eq. (|5]) (with Z — 1 of course) we have 
finally that BR((7r+ e -) -> (/i+e - )i/ M ) = 4.51 x 10 -11 . 



3.) Relativistic Effects in Muonic to Electronic Atom Decays 

We now consider (fJL~Z) — > {erZ)Vfjp e , using with obvious replacements the same 
kinematics as in Fig. 1. For simplicity we will take the nucleus to be spinless. As we 
now consider a 3-body decay, we first discuss the ranges of the relevant kinematical 
variables. Because the two constituents have extremely different masses (especially 
if Z is large), it is useful to write the masses and mg 2 of the bound states B 1 
and B 2 in the form 

m Bl = M + m 1 , m 1 = m li - E Und>1 ; 

m B2 = M + m 2 , m 2 = m e - E hin d,2 ■ (21) 

Here M is the mass of the heavy nucleus. In the rest frame of the decaying atom, 
the momentum transfer \q \ to the neutrino-pair is kinematically restricted to the 
range 

0< |g| < (m 2 Bl -m% 2 )/{2m Bl ) = {m x - m 2 ) + 0{l/M) . (22) 

For a given value of \q\, the fourth component q° as seen in the rest frame, and the 
square of the four- vector, q 2 , are fixed: 

q° = ( mi -m 2 ) +0(1/M), q 2 = [m l -m 2 ) 2 -q 2 . (23) 

Thus q 2 lies in the range [0, (mi — m2) 2 }. The momentum transfer \q \ is always 
very small compared to the masses of the bound states, so the bound states can be 
considered to be at relative rest; in this approximation our formalism will still be 
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covariant up to 0(1/M) corrections. Writing the relevant four- Fermi operator for 
muon decay in charge retention form, the matrix element for the transition reads 

4(jr 

M sr = ^m Bl m B , 2 S» , (24) 



where 



= u(p Vlt ) 7 M v(p„J (25) 

* = / S M " i] '" 4Ch - ?) r ' ,{e ' ^ - r 1 ;75 M ^ l) (26) 



= \jml + kl, t? = ^ + (*i-f) ! ■ (27) 

The muon spin r and the electron spin s are just the spins of the bound states Bi 
and B 2 . In general a Wigner rotation || must be used to relate spins in two different 
frames, but here the rotation is essentially unity because the relative velocity of B\ 
and B2 is small. 

We seek an expression for M sr that is correct to zeroth order in 1/M. This is 
found most easily by first writing the trivial identity 

S£ r = 6 aa ,S$ r ,Sr>r , (28) 

and then by rewriting the Kronecker-deltas in terms of the spinors for the bound 
states and their constituents, 

8 r t r = (2m Bl (fcJ + m (1 )) Ur'^k^u^B^O) , 



5 SS , = (2m B2 (/c°+m e )) 1 u s (B 2 ;pb 2 ) « S '(e; k x - q) . (29) 

The first relation is exact, and the second introduces errors only 0(1/M). Then Sg r 
can be written as 

^ r = (4m Bl m B2 )- 1 / 2 ^( J B 2 ;p i?2 )r^ r (S 1 ;0) , (30) 

where T M may be derived using eqs. (|26|), ( p8|) , and (|29|) . By repeatedly using the 
Dirac equations 

^°u r (B 1 ;6)=u r (B 1 ;6) , u s (5 2 ; p B2 ) 7 = « S (B 2 ; p B J + 0(1/ M) , (31) 

and dropping terms ~ g M , which vanish when contracting with the neutrino tensor 
of Eq. (^) , we arrive after some not completely straightforward algebra at 

T» = F l (q 2 )rL + F 2 (q 2 )YR + F 3 (q 2 )Y —L + F 4 (q 2 )Y—R • (32) 
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Here L = (1 — 7s)/2, and R = (1 + 7s)/2, and the form factors F(g 2 ) are given as 



with 



/M 
/ ?2 

/i.i 



d 3 A;i 
(2^)3 



^i(fci) ^2(^1 - q) 



hi 



(33) 



(*° + m M )(A; 2 + me) + q° [(1 - C)^ + m„) - C(A; 2 + m e ) 

+ (5-C)(g°) 2 -A 
(C - 5)g 2 - 
;i-C)(A;° + m /i ) + (S-C) g c 



m,, 



m. 



C(k° 2 + m e ) -(B- C)q° 
The quantities A, B, and C are 

3(£i ■ 9 f ~ q 2 kf 



_ - (fc t ■ g ) 2 

2<f 2 



2(9 



2\2 



c 



fcl ■ <f 



(34) 



(35) 



The form factors F{ may seem non-covariant, because after the integration d 3 ki the 
variables q° and \q\ remain as well as the square of the four- momentum, q 2 . But the 
form factors were derived in the rest frame of B\, so q° and \q\ are functions only 
of q 2 according to Eq. fl23D . 

In terms of the quantities introduced, the matrix element can be written in the 
suggestive form 



M,, 



~V2 



u{p v „)l^Lv(p Ve ) u s (B 2 ; pb 2 ) u r (Bi, 0) 



(36) 



The calculation of the decay rate, differential in the momentum transfer \q\, is now 
standard. Expressing the bound state masses ms 1 and m# 2 in terms of M, mi, 
and m 2 as defined in Eq. (§3), one obtains 



dr 

d\q\ 
K{\q\) 



G 2 F \q\ 
12tt 3 



-K(\q\) 



+ 2 ( mi _ m2 ) 2 ] (F 2 + F 2 ) + -\ [4(mi - m 2 ) 2 - g 2 ] (F 2 + F 4 2 ) 



mf 



-6q 2 



F^ 2 + F 3 F 4 + mi m2 (iq - F 2 ) (F 3 - F 4 ) 



mf 



m. 



(37) 



where the Fj are the form factors defined in Eq. (0). Note that mi and m 2 enter 
only through their difference; this is clear from the decomposition ([H]), which is 
invariant under the change of variables M — > M + A and m« — > rrij — A. 

The wavefunctions that enter the form factors Fi in Eq. ( |33"1) are plotted in Figs. 2 
and 3 (solid line). For comparison we have also drawn the wavefunctions for a point 
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nucleus. The finite size of the nucleus affects the muon IS wavefunction at high Z, 
as shown in Fig. 2. The finite size also affects the ultra- relativistic momentum 
tail of the electron IS wavefunction, as shown in Fig. 3. The calculation of these 
wavefunctions is discussed in detail in the appendix. Briefly, we have approximated 
the shape of a nucleus of atomic number A as a homogeneously charged sphere of 
radius r = 1.3A 1//3 fm. In position space we solve the Dirac equation inside and 
outside r ; the condition that these solutions match at r = r determines the ground 
state wavefunction. We then take its Fourier transform. 

From the kinematics (see Eq. fl22D) it is clear that the muon IS wavefunction 
is tested to momenta the order of m M . As shown in Fig. 2, relativistic effects are 
moderate. However, the finite size of the nucleus enhances the low momentum part 
of the muon wavefunction, which is fortunate because for lower momenta the electron 
wavefunction is large, and so the overlap increases. The electron IS wavefunction 
is also tested up to momenta the order of m M and so the ultra-relativistic tail of the 
Dirac wavefunction is important. As shown in Fig. 3 the finite size of the nucleus 
diminishes the tail, and so the overlap with muon wavefunction decreases. At high 
Z we find that while the shape of the spectrum in \q\ is altered, the increase and the 
decrease in the total rate due to the finite nuclear size balance remarkably, giving 
the same predicted total rate for atomic alchemy as does a calculation using a point 
nucleus. 

Using the same approximations we also calculate the free-electron decay rate 
T e j ree for (Zfi~) decays. We write the mass of the bound state B x as m Bl = 
M + 7m M , so that ^m^ is the total energy of the muon (see Eq. (p2|) . We find 

r eJree = rV (L- 1 ) , (L~ 1 ) = f |%hM£i)l 2 />2 m " = , (38) 

where T = G^m^/(1927r 3 ) is the decay width of a free muon, and (L -1 ) can be 
interpreted as the mean inverse Lorentz factor representing the slowing of the muon 
decay rate due to its orbital velocity. Numerically, we have (L^ 1 ) = 0.96, for Z = 80 
and A = 200; and (L _1 ) = 1.00, for Z = 10 and A < 20. For Z = 80 and A = 200, 
numerically 7 is 0.91; for a point nucleus with Z = 80 note that 7 and (L~ l ) are 
significantly smaller, 7 = (1 — (Za) 2 ) 1 ^ 2 = 0.81 and (L^ 1 ) = 0.85. However, for 
Z m 80, muon capture dominates the free-electron decay by a factor of about 30 
0. The branching ratio is now obtained by numerically integrating the spectrum 
in Eq. Q3"7| ) and dividing by T tot = 30 x T e j ree , with T e j ree given in Eq. Q3"8|). Using 
Dirac wavefunctions, the results for Z = 80 are 

BR[(Zfj,-) -> (Ze-)u^ e ] = 1.19 x 10" 6 , (finite nucleus, A = 200) , 
BR[{Z^-) -> (Ze-)v^ e ] = 1.19 x 10~ 6 , (point nucleus) . (39) 

We note that the effects of the finite size of the nucleus do influence the shape of 
the momentum spectrum dT/d\q\, as shown in Fig. 4, but that they leave the total 
rate essentially unchanged even for Z as large as 80. As the finite size effects are 
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expected to be smaller for smaller Z, we can calculate the rate for Z = 10 using the 
Dirac wavefunctions for a point nucleus. We get now, assuming « equal rates for 
capture and free-electron decay 

BR[(Zfi~) -> (Ze~)z/ M z7 e ] = 1.25 x 10" 9 , (point nucleus, Z = 10) . (40) 

To appreciate the relativistic enhancement due to the Dirac wavefunctions, we 
also show the result obtained with the Schrodinger wavefunctions for a point nucleus: 

BR[{ZiT) -> (Ze~)z/ M z> e ] = 2.32 x 10" 8 , (Z = 80) , 

BR[(Zfi~) -> (Ze _ )i/ M z/ e ] = 9.32 x 10~ 10 , (Z = 10) . (41) 

The relativistic enhancement for Z = 80 is a remarkable factor of 50. 

In Fig. 4 the decay distribution e£T/<2|g|, divided by the total width T tot , is given 
for Z = 80. For illustration we also give the spectra predicted using the Schrodinger 
and Dirac wavefunctions for a point nucleus. Comparing these two curves shows that 
relativistic effects not only change the overall normalization but cause the spectrum 
to peak at higher momenta. The finite nuclear size causes the shape of the spectrum 
as calculated for a point nucleus to narrow. This can easily be understood by looking 
at the muonic wavefunction in Fig. 2. In Fig. 5 the same distribution is shown for 
Z = 10; as might be expected both relativistic effects and the effect of a finite nuclear 
size are small. Transitions of the form (Zfi) — > (Ze) + vv have as their signature 
a bound state recoiling with a large momentum the order of m M . The momentum 
distribution for heavy atoms ( Z = 80 ) peaks at about 20 MeV as seen in Fig. 4. 

4.) Conclusions 

In this paper we have shown that transitions between electromagnetic bound 
states can be calculated reliably using relativistic atomic wavefunctions. The rates 
and also the shapes of the spectra depend drastically on relativistic corrections; 
for heavy muonic atoms these enhance the branching ratio to a sizeable value of 
~ 10~ 6 from ~ 10 -8 (see Eqs. ( |3~9~D and (f4T|)). The QED analysis of these alchemy 
transitions illustrates some of the physics of the relativistic wavefunctions that must 
invariably enter the QCD analysis of the corresponding exclusive electroweak decays 
of hadrons. 

Prospects for an experimental test of atomic alchemy are dim but not hopeless. 
While the branching ratio appears to be too small to be detectable for (ire) or for 
light muonic atoms, it may be possible to detect it for heavy muonic atoms. One 
promising approach is to inject and capture ~ 5 kev muons in a cyclotron trap Q 
containing ~ 0.1 bar of a noble gas, which we will here assume to be neon [Z = 10). 
The muons will come to rest Q in ~ 1 /is and be captured by a neon atom. The 
muons will rapidly cascade to the IS* state, mostly by the ejection of atomic electrons. 
In neon the mean number of electrons in the K shell has been measured [EJ to lie 
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between 0.07 and 0.68 by the time the muon has fallen to states with principal 
quantum number n = 5, in agreement with simple estimates || of 0.25, and will 
certainly diminish further as the muon continues to the IS" state. Total ionization 
of an atom has been observed in the cascades of antiprotons in neon, argon, and 
krypton 0. The electron K shell will therefore be vacant to receive the electron 
from atomic alchemy. Only 0.31 of the muons in the IS" state of neon will be lost 
to muon capture |§, so most of the captured muons will be available to decay by 
atomic alchemy. At a pressure of ~ 0.1 bar, using estimates in |J for the velocity 
of muonic neon ion (~ 10 5 cm/s) and for the cross section for electron capture for 
Ne 9+ on Ne (~ 3 ■ 10~ 15 cm 2 ), the time taken for a bare neon atom to capture an 
electron is the order of 1 /is, so most of the captured muons will decay before the 
electron K shell can refill and block the atomic alchemy. It has been suggested [Q 
that cyclotron traps, admittedly with higher gas pressures, can be built capable of 
stopping ~ 10 7 fi"/s, so a branching ratio as small as ~ 10~ 9 might be accessible. 

Even if such a system could be realized the experimental difficulties would be 
daunting. The signature for atomic alchemy is the appearance of 20 Ne 9+ ions with a 
distribution dN/dq 2 extending to a momentum the order of m M . These ions will be 
difficult to distinguish from 20 F 9+ ions with a momentum of m M made by internal 
conversion. Nor will it be easy to extract the ions for analysis because in 0.1 bar of 
Neon the distance a bare ion can fly before capturing an electron is only ~ 1 mm, and 
the gas pressure cannot be lowered without letting the muons in the trap decay in 
flight before they can be captured. However, measurement of the spectrum dN / dq 2 
would be valuable; it directly reflects the local structure of the atomic momentum 
space wavefunction and the beta decay spectrum. This appears to be one of the few 
ways in which one can directly study the relativistic tail of bound state wavefunctions 
in QED. 

For induced reactions, such as 7 (d + e~) — > (p + e~)n, one will observe a nearly 
monoenergetic final state hydrogen atom recoiling against an outgoing neutron. In 
principle one can use the doppler-shifted radiation from the outgoing atoms as a 
precise calibration of the nuclear scattering kinematics. One could also analyze the 
spin states of the outgoing atomic system and its hyperfine spectrum as a probe of 
the spin state of the final state nucleus. 

We would like to thank P. Truol, R. Engfer, LB. Khriplovich, and M. Strikman 
for stimulating discussions. 
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Figure captions 
Figure 1 

General diagram for a weak alchemy transition B\ — > B 2 + X as discussed in the 
text. The momenta of the particles are indicated in parentheses. 

Figure 2 

Schrodinger and Dirac wavefunctions (multiplied by \k\) for a {Zfi~) atom with 
Z=80. 
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Figure 3 

Schrodinger and Dirac wavefunctions (multiplied by \k\) for a (Ze~) atom with 
Z=80. 

Figure 4 

Decay spectrum dT/d\q\ divided by T to t = 30 x Y e j ree (Eq. (|38|)) for Z=80 (see text). 
The full line is obtained by using Dirac wavefunctions which take into account the 
finite size of the nucleus. The dashed (dashed-dotted) line corresponds to point-like 
Dirac (Schrodinger) wavefunctions. The curve based on Schrodinger wavefunctions 
is multiplied by a factor 50. 

Figure 5 

Decay spectrum dT/d\q \ divided by T tot = 2 x T e j ree (Eq. (RST)) for Z=10 (see text). 
The full (dashed) line corresponds to point-like Dirac (Schrodinger) wavefunctions. 

Appendix: Dirac wavefunctions 

In atomic alchemy (Zfi) — > (Zejv^v^ the electron wavefunction is probed at mo- 
menta the order of m^. Schrodinger wave functions are not appropriate and the 
Dirac wavefunctions || must be used. Recoil corrections due to the finite nuclear 
mass are extremely small, so the Dirac wavefunctions for an electron in the field of 
an infinitely heavy nucleus will describe adequately the high momentum tail of the 
electron (and muon) wavefunction. At large Z the Bohr radius of the muon IS state 
is less than the nuclear radius, so the effect of the finite size of the nucleus on the 
muon wavefunction must obviously be included; the effect of the finite size is also 
important on the high-momentum tail of the electron IS wavefunction. 

We sketch how the momentum-space wavefunctions ipi{k) and ipzik) that appear 
in the form factors in Eq. (^) are extracted from the usual 4-component Dirac 
wavefunctions in position space. Because everything can be worked out analytically 
for a point nucleus, we describe this case first. 

Point nucleus 

From the literature (e.g. from || on page 55) we take the (four component) ground 
state wave function in position space, ipn=i,j=i/2,j z =i/2( r , ^ 4>)- F° r brevity we write 
this as $(af). Taking the Fourier transform, we get the wavefunction in momentum 

space, 

$(£) = / d 3 x$(x) exp(-ik ■ x) , (42) 
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f(k) 



a • k 



1 




— |^| •> 



(43) 



with 



N 



(sin p + ak cos p) 



(44) 



A;(l + a 2 A; 2 ) 1+7/2 

,2n m(1 2 72^+7/2 K 1 + ( X + ^ a2e ) sin /> - ^ cos P } ( 45 ) 
7«r (1 + cr/r) " 



1 



mZa 



, 7 



1 - (Za) 2 , p = 7 atan(aA;) , iV = 2 7+1 T(l + 7) 



07r(l 



\ r(i + 27) 



(46) 

Here m denotes the reduced mass of the system, which we take to be identical to 
that of the muon or the electron because we work to lowest order in 1/M. The 
corresponding energy eigenvalue is then E = m^. For Z = 80 the numerical value 
of 7 is 0.81. 

The Dirac wavefunction for a bound state of course projects onto plane waves 
of both positive and negative energies; the latter waves correspond to antiparticles. 
We therefore define the relevant wavefunction (to be used in calculating the form 
factors) by the projection on positive energy plane waves. Thus we expand $(A;) in 
terms of spinors u r (k) and v r (—k), writing 



<b(k) = £ 

r 

If j z = 1/2 we get 



A r ik) u ^ +B :x-k) v ^' k 



l +l/2 



k° + m 
2k° 



f(k) 



2k° 



k 



k° 



k 2 + m 2 



(47) 



B 



+1/2 



k° + m 

t k° + m (k 1 + ik 2 ) ( k 
2k° k \k° + m 

'k° + m k 3 ( k 



g(k)j ; A_ 1/2 = 0; 
f(k)-g(k)) 



(48) 



2k° k \k° + m 



f(k)-g(k) 



Here A r (k) is the probability amplitude to find an electron with momentum k and 
spin r in the atom, while B*(—k) is the probability amplitude to find a positron 

— * 

with momentum k and spin r; the latter amplitude arises from the creation of e + e _ 
pairs on the nucleus. Because the wave function in position space is normalized 
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as / d 3 x $ + (x ) $ (x ) = 1, the Fourier transform is automatically normalized as 
/ d 3 k/(2ir) 3 & + (k) $(&) = 1. Therefore A r (/c) and B*(—k) are normalized so that 

/ (^E{l^)| 2 + I^)| 2 } = l • (49) 

The integral / d 3 k/(2ii) 3 J2 r \Br{k)\ 2 gives the probability to find a three particle 
Fock state (e + e~e~) in the atom. Even for Z = 80 this fraction is tiny (« 0.2%), 
so we only consider the one- Fock contribution characterized by A r (k). We mention 
that if we consider the atom with j z = —1/2, we get A + y 2 (k) = 0, and A^i/ 2 {k) is 
identical to A + i/ 2 (k), given in Eq. ( pE8"D for j z = +1/2. The wavefunction denoted 
in the text as ijj(k) is therefore given in the relativistic case by 

^j(k) = A +1/2 (k) . (50) 

Because the coefficients B r are very small, the effects of antiparticle production are 
small, and we get essentially the same form factor as we would have got had we 
naively computed the simple overlap of the wavefunctions in position space. 



Finite-size nucleus 



As usual (see e.g. flip]), we model a nucleus of atomic number A as a homogeneously 
charged sphere of radius r = 1.3 A 1 / 3 fm; numerically r ~ 7fm for A = 200. Inside 
the sphere the potential has the form V(r) = —(Za/r Q ) ■ (3 — r 2 /tq)/2, and outside 
we have V(r) = —Za/r. In the notation of Landau-Lifschitz [JTTJ| we write the 
four-component Dirac function $(f) as 



$(f) = N 



\ 




where the constant N is chosen such that / d 3 r <3> + (r) $(r 
tions for / and g are 

f'(r) - (E + m-V(r))g(r) = 

2 

g'(r) + -g(r) + (E - m- V{r)) f(r) = 



(51) 



1. The radial equa- 








(52) 



which are solved separately in the two regions r < ro and r > ro for an arbitrary 
constant E < m. In the outer region r > tq there is (up to an overall constant) 
exactly one solution (/, g) that is square integrable at r = oo. It is given by (see 
Landau-Lifschitz [ |TT|] ) 



f(r) 



Vm + E exp(-Ar) r^ 1 {Q x + Q 2 ) 
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g(r) = —y/m — E exp(— Ar) r 7 1 (Qi — Q2) 

ry 771 

Qi = C/( 7 --?-,2 7 + l,2Ar) 

„ , Zam*. , ZaE , , 

Q 2 = )E% + 1-— ,2 7 +l,2Ar) 

A = Vm 2 - £ 2 , 7 = ^1 - (Za) 2 . (53) 

The hypergeometric function U is denned and discussed in detail in chapter 13 of 
Abramowitz-Stegun |T2]; it is also related to Whittaker functions which are numer- 



ically accessible in the CERN-library. 

In the inner region (r < r ) a simplified solution for the (massless) electron 
has been obtained by Khriplovich ||13|| : however, the substantial mass of the muon 
requires a more complete treatment. It turns out that the Dirac equation has (again 
up to an overall constant) exactly one solution (/, g) that is square integrable at 
r = ro; in the present case it has a Taylor series expansion around r = 0. While / 
starts as a constant, g begins with a term linear in r. The coefficients of the power 
series expansions of / and g may be defined recursively. 

The inner and outer solutions must satisfy the matching condition 

{f/g)r^r - = (f/g)r-+r 0+ (54) 

in order to be solutions of the complete equation. This can only hold for certain 
values of E; these are just the eigenvalues. For Z = 80 and A = 200 the lowest 
eigenvalues are E = 0.908m M and E = 0.811m e for the IS* states of the muon and 
the electron, respectively. 

The Fourier transform $(&) is then defined and written in terms of f(k) and 
g{k) precisely as in eqs. (^) and (f|3|). Proceeding through the same steps as for 
a point nucleus, one finally gets the wavefunction ip{k) (compare with Eq. (|50|) ) in 
a numerical form. This wavefunction shown for the [Z[i) and the [Ze) atoms in 
Figs. 2 and 3, respectively. 
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